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del género Boloria Moore, 1900 en la Península Ibérica 

(Lepidoptera: Nymphalidae) 
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Resumen 


El cambio climático está afectando a los rangos de distribución de los Lepidoptera, lo que podría afectar más a 
especies propias de montaña en zonas templadas como el sur de Europa. El presente estudio pretende analizar el gé- 
nero Boloria Moore, 1900 en la Península Ibérica, ya que comprende varias especies de alta montaña, mediante el 
uso de modelos predictivos bajo diferentes escenarios climáticos (para los años 2020, 2050 y 2080). Los resultados 
obtenidos fueron significativos (p=0,05), con valores de AUC (Area Under a Receiver Operating Characteristic 
(ROC) Curve) superiores a 0,9 lo que indica que las predicciones fueron mejores que las obtenidas al azar. Las va- 
riables climáticas que mostraron mayor influencia en los modelos fueron la precipitación del trimestre más cálido y 
la temperatura media del trimestre más seco. Nuestros modelos muestran una tendencia general a la reducción del 
área de distribución favorable en todos los escenarios para todas las especies, llegando éstas incluso a desaparecer en 
2080 en algunos casos, como B. pales y B. napaea. El cambio climático puede estar afectando negativamente a este 
género, por lo que sería pertinente hacer un seguimiento de estas especies particularmente sensibles para comprobar 
en un plazo más corto si estas tendencias se confirman. 

PALABRAS CLAVE: Lepidoptera, Nymphalidae, Boloria, cambio climático, AUC, MaxEnt, precipitación, 
temperatura, Península Ibérica. 


Predicting climate change impacts on Lepidoptera species distribution. 
The case of genus Boloria Moore, 1900 in the Iberian Peninsula 
(Lepidoptera: Nymphalidae) 


Abstract 


Climate change is affecting the geographic ranges of Lepidoptera. Such changes are presumably more marked 
in mountain species in areas such as southern Europe. The present study tries to analyze the genus Boloria Moore, 
1900, in the Iberian Peninsula since it comprises several mountain species. Different predictive models were applied 
considering all the Iberian species in the genus and three different climatic scenarios (the years 2020, 2050 and 
2080). The results obtained were significant (p= 0.05) and AUC values (Area Under a Receiver Operating 
Characteristic (ROC) Curve) were higher than 0.9. AUC values obtained indicate that the predictions were better 
than random predictions. The most significant variables in the models were “raifall of the warmest quarter of the 
year” and “mean temperature of the driest quarter”. Our models show that there is a general tendency towards a 
decrease of favorable distribution areas in all future scenarios for all species, going as far as extinction in 2080 in 
some cases, such as B. pales and B. napaea. Climate change may be affecting this genus negatively so that it would 
be appropriate to closely monitor these particularly sensitive species to determine, within a shorter time frame, 
whether these trends are confirmed. 

KEY WORDS: Lepidoptera, Nymphalidae, Boloria, climatic change, AUC, MaxEnt, precipitation, temperature, 
Iberian Peninsula. 
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Introducción 


Lepidoptera es uno de los órdenes de insectos con mayor número de especies descritas (GU- 
LLAN & CRANSTON, 2005). En los últimos años se ha observado una retracción del área ocupada 
por ciertas especies (WILSON et al., 2007). Estas disminuciones pueden estar causadas por cambios 
en el uso del suelo y el manejo del territorio (NEW et al., 1995), por la introducción de especies 
alóctonas (DIDHAM et al., 2007) o por el cambio climático (PARMESAN et al., 1999; HILL et al., 
2002; THOMAS et al., 2004) entre otros. Todo ello se traduce en efectos tales como aislamiento, 
fragmentación y pérdida de hábitats favorables, ya documentados para diversas especies del grupo 
(VAN SWAAY & WARREN, 1999; WARREN et al., 2001; STEFANESCU et al., 2011). 

De los factores anteriormente mencionados, el efecto del calentamiento global en especies de 
diferentes grupos taxonómicos es actualmente objeto de especial atención (por ejemplo THOMAS et 
al., 2004; ARAÚJO et al., 2011; FELICÍSIMO, 2011; STEFANESCU et al., 2011). En las últimas 
décadas ha tenido lugar una alteración importante de las condiciones atmosféricas que no parecen re- 
lacionarse con la “variabilidad climática natural” (cambios en el clima que se han ido dando en la 
historia de la tierra por variaciones de ciclos dinámicos naturales). De ello puede deducirse que el 
cambio climático es debido a otras causas además de las naturales, siendo una de las principales la 
actividad humana, aspecto sobre el que existe ya un amplio consenso (ORESKES, 2004) por ser el 
hombre el responsable del aumento de emisión de gases de efecto invernadero (IPCC, 2001) que tie- 
ne consecuencias en la circulación del aire atmosférico, en el patrón de precipitaciones e incluso en 
el nivel del mar por deshielo de los glaciares, estableciéndose una destacada tendencia de aumento 
de la temperatura global. Si esta tendencia permanece en el tiempo, podría convertirse en una seria 
amenaza para las poblaciones de Lepidoptera en Europa (VAN SWAAY & WARREN, 1999; THO- 
MAS et al., 2004). En comparación con muchos vertebrados, las poblaciones de estos insectos pare- 
cen especialmente proclives a oscilaciones anuales provocadas por factores demográficos relaciona- 
dos con variables medioambientales (POLLARD, 1988; ROY et al., 2001). Ciertas características, 
como el hecho de poseer un ciclo de vida corto, confieren a los Lepidoptera gran sensibilidad ante 
variaciones en su hábitat; ello los convierte en buenos indicadores de cambio climático (POLLARD, 
1988; PARMESAN, 2003; STEFANESCU ef al., 2003). Por ello el análisis de los factores que condi- 
cionan su distribución geográfica constituye una herramienta útil para comprender las tendencias de 
la biodiversidad a gran escala (THOMAS, 2005; VAN SWAAY et al., 2006). 

Las distribuciones de las mariposas diurnas están condicionadas por multitud de factores, ya se- 
an climáticos (HAWKINS et al., 2003) como los que se trata de comprobar su efecto en este estudio, 
o de otro tipo como ambientales o bióticos (p.ej. competencia y depredación). Debido a esta comple- 
jidad, y por la dificultad de delimitar las áreas de distribución de las especies con suficiente preci- 
sión, en especial en zonas insuficientemente estudiadas, se hace tentador recurrir a modelos predicti- 
vos para estimar las distribuciones potenciales (ROMO et al., 2006; JIMENEZ-VALVERDE et al., 
2008). Estos modelos se utilizan porque indican la idoneidad del habitat para el desarrollo de pobla- 
ciones de una especie concreta a partir de sus datos de presencia (o ausencia) (ARAÚJO & GUISAN, 
2006; FERRIER & GUISAN, 2006); permiten, además, evaluar la influencia de diversos factores cli- 
máticos sobre la distribución de las especies. Entre los modelos y escenarios aplicables, los más co- 
munes son las proyecciones de las distribuciones a los años 2020, 2050 y 2080, eficaces para deter- 
minar distribuciones potenciales admitiendo cierto margen de error en los datos (PEARSON & 
DAWSON, 2003; BENITO DE PANDO & PEÑAS DE GILES, 2007; PINEDA & LOBO, 2009). La 
validación es en cualquier caso problemática, ya que no se puede comprobar fácilmente algo que to- 
davía no ha pasado (ARAÚJO et al., 2005). En muchos casos tendremos que contentarnos con una 
valoración independiente de cualquier conjunto de datos. Como los Lepidópteros de ambientes mon- 
tañosos parecen ser más susceptibles al cambio climático (BENISTON et al., 1997), en este estudio 
se seleccionó el género Boloria Moore 1990 (superfamilia Papilionoidea, familia Nymphalidae), re- 
presentado en la Península Ibérica por seis especies que son características de altitudes y latitudes 
elevadas. Para la realización de los modelos se utilizó el programa informático MaxEnt, que presenta 


268 SHILAP Revta. lepid., 41 (162), junio 2013 


PREDICCION DE LOS IMPACTOS DEL CAMBIO CLIMATICO EN LA DISTRIBUCION DE LEPIDOPTERA DEL GENERO BOLORIA 


ventajas cuando sólo se poseen presencias como datos de partida (PHILIPS et al., 2006) y debido a 
que se encuentra entre los métodos más fiables de modelización (ELITH et al., 2006; BENITO DE 
PANDO & PEÑAS DE GILES, 2007) tanto para predicciones basadas en condiciones actuales como 
futuras (HIJMANS & GRAHAM, 2006). 

Los objetivos de este estudio consisten en: 1) elaborar mapas de distribución potencial de las es- 
pecies pertenecientes al género Boloria a partir de datos actualizados de distribución conocida; y 2) 
determinar la posible influencia del cambio climático en la distribución de este género en el futuro 
mediante dos escenarios climáticos diferentes utilizando las proyecciones generadas para los años 
2020, 2050 y 2080. 


Material y métodos 
ÁREA DE ESTUDIO, REPRESENTACIÓN GEOGRÁFICA Y GRUPO TAXONÓMICO 


El área de estudio lo constituye la Península Ibérica, utilizándose para la representación de la 
distribución de las especies el sistema de proyección de coordenadas UTM (Universal Transverse 
Mercator) con redes de malla de 10 km de lado. 

El género Boloria (Papilionoidea, Nymphalidae) está representado por seis especies peninsula- 
res: Boloria dia (Linnaeus, 1767), B. eunomia (Esper, 1800) B. selene ([Denis & Schiffermiiller], 
1775), B. euphrosyne (Linnaeus, 1758), B. pales ([Denis & Schiffermiiller], 1775) y B. napaea 
(Hoffmansegg, 1804). Para la nomenclatura de las especies se utilizó el catálogo de VIVES-MORE- 
NO (1994). Todas ellas se distribuyen por la región septentrional de la Península Ibérica, alcanzando 
en ocasiones el centro de la misma (GARCÍA-BARROS et al., 2004). Presentan un rango altitudinal 
de entre 900 y 2.900 m, aunque algunas especies como B. euphrosyne o B. dia pueden encontrarse 
también a altitudes más bajas. B. eunomia, B. pales y B. napaea son las especies que alcanzan rangos 
altitudinales más elevados. 

Los datos de distribución proceden de la base de datos de presencias de mariposas compilados 
por GARCÍA-BARROS et al. (2004). Estos datos se actualizaron para el género Boloria con diversas 
publicaciones hasta el año 2010. La base de datos utilizada cuenta con 4.711 registros, de los cuales 
1.979 eran referentes a la especie B. dia, 1.242 a B. selene, 985 a B. euphrosyne, 56 a B. eunomia, 
431 a B. pales, y 18 a B. napaea. 


VARIABLES CLIMÁTICAS 


Para la realización de los modelos de distribución potencial actual se utilizaron las 19 variables 
bioclimáticas (BIOCLIM) del presente procedentes de la base de datos de WorldClim -disponible en 
http://www.worldclim.com y descritas por HIJMANS et al. (2005)-, con una resolución de 5 minutos 
(que es la resolución que más se asemeja a la de los datos de presencia de las especies medidas en 
cuadrículas de 10x10 km de lado). Estas variables bioclimáticas son el resultado de recopilaciones de 
medidas mensuales de distintas variables de temperatura y precipitación llevadas a cabo entre el año 
1950 y 2000 por un gran número de estaciones meteorológicas repartidas por todo el mundo, con el 
objetivo de conseguir unas capas ambientales fiables y biológicamente más significativas (HIJ- 
MANS et al., 2005). Representan tendencias anuales, estacionalidad y factores ambientales extremos 
o limitantes. 

Para la realización de los modelos de distribución potencial influidos por el cambio climático, 
se obtuvieron con una resolución de 5 minutos las mismas 19 variables bioclimáticas proyectadas al 
futuro del CIAT-GCM (Portal de datos del Centro Internacional de Agricultura Tropical, http://ccafs- 
climate.org/download_sres.html). Se eligieron las variables climáticas generadas por el Modelo de 
Circulación General (MCG) de la atmósfera HCCPR HadCM3 (Hadley Centre Coupled Model, ver- 
sion 3) desarrollado en el Hadley Centre for Climate Prediction and Research (HCCPR, Reino Uni- 
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do). Los datos fueron calculados por el HCCPR a partir de representaciones relativas al período 
1961-1990 y distribuidos por el IPCC (Intergovernmental Panel on Climate Change). 

Debido a que cada proyección de cambio climático depende de las consecuencias de la activi- 
dad humana que se realizarán a lo largo del tiempo, y no se puede asegurar como acertada una única 
proyección para determinado año, existen diferentes escenarios denominados SRES (Special Report 
on Emission Scenarios) para evaluar el cambio climático y realizar los modelos. Los SRES seleccio- 
nados en este caso fueron el A2 y el B2. El A2 representa una proyección de emisiones de CO, co- 
rrespondientes a un elevado consumo energético (situación más pesimista en cuanto a gases de efecto 
invernadero se refiere) y prevé un cambio de temperatura en la superficie del aire entre los años 2071 
y 2100 de 3,0°C como promedio (en comparación con el período de 1961 a 1990; IPCC, 2001). El 
B2 representa una proyección de emisiones correspondientes a un bajo consumo energético y prevé 
un calentamiento más leve, con un cambio de temperatura entre 2071 y 2100 de 2,2*C (IPCC, 2001). 
Ambos escenarios proyectan un aumento de los promedios mundiales de vapor de agua, evaporación 
de la misma y precipitaciones, aunque a escala regional, concretamente en el área de estudio ibérico 
que nos concierne, se observan aumentos en latitudes elevadas, y disminuciones leves en latitudes 
medias. En resumen, estos escenarios son una gama de imágenes de posibles situaciones futuras, que 
se construyen a partir de diversos factores (como por ejemplo contaminación de gases, uso del suelo 
o consumo de combustibles fósiles) y del modelo HadCM3, con los que puede simularse el clima 
hasta el año 2020, 2050 ó 2080. 


ANÁLISIS PREVIOS 


Las variables BIOCLIM podían mostrar un alto nivel de correlación entre ellas que distorsiona- 
ría los resultados obtenidos y aumentaría la redundancia de datos, por lo que se llevó a cabo una se- 
lección previa de las variables a utilizar. Se extrajeron los valores de estas capas climáticas con el 
programa informático Arc Gis versión 9.2 (ESRI, 2006), y a partir de ellos se calculó el nivel de co- 
rrelación entre pares de variables en el área de estudio analizada gracias al programa estadístico 
SPSS versión 15.0 (SPSS, 2006). Siguiendo a RISSLER & APODACA (2007) cuando dos variables 
compartían un coeficiente de correlación de Pearson mayor que 0,8, se seleccionó la variable más 
significativa a nivel biológico para estas especies o aquella que resultara de más fácil interpretación. 
Como resultado, las variables climáticas no correlacionadas entre sí con las que se realizaron los mo- 
delos predictivos de este trabajo fueron: Biol (temperatura media anual), Bio2 (rango de temperatura 
medio mensual, obtenido como la diferencia entre la máxima y la mínima temperatura de cada mes), 
Bio3 (Isotermalidad), Bio8 (Temperatura media del trimestre más húmedo), Bio9 (Temperatura me- 
dia del trimestre más seco), Bio12 (Precipitación anual) y Bio18 (Precipitación del trimestre más cá- 
lido). 


MODELIZACIÓN: MAPAS DE DISTRIBUCIONES POTENCIALES 


Los modelos de distribución potencial se elaboraron con el programa Maxent versión 3.3.3 
(http://www.cs.princeton.edu/-schapire/maxent/). Este aplica el principio de máxima entropía para 
calcular la distribución geográfica más probable para una especie. De este modo, el modelo de pre- 
dicción de cada especie se materializa en un mapa que representa las áreas adecuadas o no para la 
presencia de dicha especie. Una explicación detallada de cómo funciona el programa se encuentra en 
PHILLIPS er al. (2006) y PHILLIPS & DUDIK (2008). 

Como formato de salida de los resultados se eligió el formato logístico por ser de más fácil in- 
terpretación (PHILLIPS & DUDIK, 2008), ya que pueden ser interpretados como probabilidad de 
presencia de la especie condicionada a las características ambientales. Se consideraron los paráme- 
tros por defecto (máximo número de iteraciones: 500, límite de convergencia: 0,00001, y 10.000 
puntos de fondo) como recomiendan PHILLIPS et al. (2006). Se realizaron 100 réplicas por modelo, 
utilizándose la media de cada uno de ellos para la visualización de los resultados, y el umbral de má- 
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xima probabilidad de presencia se ajustó al que representaba mejor los datos de presencia conocida 
de la especie. En todos los casos, un 25% de los registros de la muestra fueron utilizados para medir 
la calidad del modelo, y el 75% restante para realizar el modelo, tal y como sugieren diversos autores 
(PAWAR et al., 2007; ALBA-SÁNCHEZ et al., 2010). Además de esta práctica de validación sim- 
ple, se aplicó el método de validación cruzada a las diferentes réplicas de los modelos 

Con esta metodología se generaron para cada especie un modelo predictivo de la distribución 
potencial para el presente y tres proyecciones para los años 2020, 2050 y 2080 para cada uno de los 
dos escenarios de cambio climático (A2 y B2). 

Con el programa DIVA-GIS (HIJMANS et al., 2004) se calculó el número de cuadrículas predi- 
chas por los modelos por encima del umbral de máxima probabilidad, para poder valorar de manera 
cuantitativa los cambios acaecidos en las distribuciones de las especies. 


EVALUACIÓN DE LOS MODELOS PREDICTIVOS 


Para evaluar o medir la bondad de los modelos creados, se utilizó el valor AUC (Area Under a 
Receiver Operating Characteristic (ROC) Curve). Este valor indica la capacidad de discriminación 
del modelo (PEARCE 4 FERRIER, 2000). El estadístico AUC puede tomar valores entre 0,5 y 1. 
Un valor de 0,5 indicaría que el modelo obtenido no es mejor que uno realizado al azar, un AUC de 
0,5 a 0,7 corresponde a un modelo de baja precisión o escasa discriminación, mientras que un valor 
mayor de 0,7 corresponde con un modelo de elevada precisión o alta discriminación (PEARCE & 
FERRIER, 2000; NEWBOLD et al., 2009). 

Además se calculó la significación estadística de la predicción usando una prueba binomial de 
omisión para cada especie consistente en once tests binomiales como complemento para la evalua- 
ción del modelo: valor fijo acumulado 1; valor fijo acumulado 5; valor fijo acumulado 10; el mínimo 
de puntos de presencia de entrenamiento; el percentil 10 de los puntos de presencia de entrenamien- 
to; prueba de igualdad de la sensibilidad y especificidad de entrenamiento; valor máximo de la sensi- 
bilidad más la especificidad de entrenamiento; prueba de igualdad de la sensibilidad y especificidad 
de prueba; valor máximo de la sensibilidad más la especificidad de prueba; balance de los puntos de 
omisión de entrenamiento, área predicha y valor umbral; comparación de la entropía de las distribu- 
ciones originales y consideradas a un umbral determinado (PHILLIPS et al., 2006). Para que la pre- 
dicción se considere admisible se requiere que todos ellos sean significativos con valores de p<0,01. 

El valor del umbral logístico de la prueba de igualdad de la sensibilidad y especificidad de los 
valores de validación se utilizó para calcular porcentajes de clasificación que ayudaran a corroborar 
los resultados obtenidos por el estadístico AUC. Estos porcentajes se calcularon como el número de 
localidades de prueba con valores logísticos de probabilidad mayores que este valor umbral dividido 
por el número total de localidades de prueba (BALDWIN & BENDER, 2008). 


Resultados 
ÁREAS DE DISTRIBUCIÓN POTENCIAL 


1.1. Distribución potencial actual de las especies del género Boloria 

Las predicciones realizadas se ajustan bastante al rango de distribución actual de las especies 
(Figura 1). Se sitúan sobre las principales cordilleras de la zona norte de la Península Ibérica alcan- 
zando, en determinados casos como B. dia o B. selene cordilleras más centrales. 

A pesar de que en muchos casos hay un patrón de distribución coincidente situado en la cordi- 
llera Cantábrica y los Pirineos, en el caso de B. napaea la predicción se restringe solamente a la zona 
pirenaica donde se sitúan las poblaciones conocidas. 

Las variables climáticas que contribuyeron en mayor medida a la elaboración de los modelos 
fueron la precipitación del trimestre más cálido (Bio18) y la temperatura media del trimestre más se- 
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co (Bio9), seguidas por la temperatura media anual (Biol) y en el caso de B. selene la precipitación 
anual (Bio12) (Tabla D. 


Especie Variable bioclimática Porcentaje de 
de mayor contribución contribución (%) 
B. dia Bio 18 52,4 

Bio 9 20,0 

B. eunomia Bio 1 55,8 

Bio 18 21,6 











B. selene Bio 9 53,1 


Bio 12 15,6 
B. euphrosyne Bio 1 33,7 
Bio 9 29,0 
B. pales Bio 1 67,7 
Bio 18 24,7 


B. napaea Bio 18 83,4 
Bio 9 7,3 




















Tabla I.- Porcentaje de la contribución de las dos variables bioclimáticas más utilizadas por el programa MaxEnt en 
la elaboración de los modelos predictivos de cada especie de mariposa del género Boloria. Biol (temperatura media 
anual), Bio9 (Temperatura media del trimestre más seco), Bio12 (Precipitación anual) y Bio18 (Precipitación del 
trimestre más cálido). 


1.2. Distribución potencial futura de las especies del género Boloria para los años 2020, 2050 y 
2080 

Para cada especie se obtuvieron seis modelos predictivos de la distribución potencial: uno para el 
año 2020, otro para 2050 y otro para 2080 bajo los escenarios de cambio climático A2 y B2 respectiva- 
mente (Figuras 2-7). Los umbrales elegidos como probabilidad de presencia difirieron en cada especie, 
eligiéndose en cada caso aquel que proporcionaba la predicción más parecida a la distribución conocida 
de cada una de ellas. Se consideró un umbral de 0,55 para las especies B. dia y B. selene; de 0,6 para B. 
euphrosyne, B. eunomia y B. pales; y de 0,75 para B. napaea. 

Se aprecia una pérdida generalizada del área de ocupación sobre la distribución potencial predicha 
en los modelos obtenidos a lo largo de los años en todas las especies estudiadas (Figuras 2-7; Tabla ID. 
Este patrón fue más acusado cuando el modelo se proyectó sobre el escenario de cambio climático A2, 
salvo en el caso de B. dia, B. pales y B. napaea. 


Distribución Predicción 
conocida presente 
B. dia 547 811 


B. eunomia 122 


Especie 











B. selene 





B. euphrosyne 





B. pales 





























B. napaea 





Tabla II.— Número de cuadrículas UTM que conforman la distribución conocida de cada especie del género de 
mariposas Boloria y el número de cuadrículas que se obtienen en los modelos de predicción, tanto en el presente, 
como en el futuro para los años 2020, 2050 y 2080 bajo los escenarios climáticos A2 y B2 (ver texto). Se señala el 
umbral elegido para el recuento de las cuadrículas de las predicciones según la especie. 
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En el año 2020 las diferencias con las predicciones actuales no son notables. En el año 2050 hay 
un claro patrón de concentración de posibles áreas favorables en las zonas de mayor altitud de la región 
septentrional de la península, como son los Pirineos y algunos puntos de la cordillera Cantábrica. En el 
año 2080 hay una reducción drástica de las áreas favorables, llegando a desaparecer casi por completo 
para algunas especies, como por ejemplo B. pales y B. napaea. Otras especies como B. dia, B. selene y 
B. euphrosyne, todavía presentaron áreas favorables en los Pirineos y algunas cuadrículas de la cordille- 
ra Cantábrica en el año 2080, mientras que en B. eunomia desaparecen las poblaciones de la cordillera 
Cantábrica. 

B. dia es la especie que presenta mayor cantidad de cuadrículas favorables en el año 2080 en el es- 
cenario A2, llegando a presentar hasta 204 cuadrículas y B. selene es la que presenta mayor número de 
cuadrículas en el escenario B2, con un total de 178 cuadrículas. 


EVALUACIÓN DE LOS MODELOS PREDICTIVOS 


1.1. Fiabilidad del modelo: Curva ROC y estadístico AUC 

El valor AUC obtenido en los modelos predictivos realizados para las seis especies del género Bo- 
loria fue superior a 0,9 (Tabla III), de modo que desde este punto de vista los modelos obtenidos se 
pueden considerar fiables (aceptables) y más precisos que un modelo obtenido al azar. 


Especie Umbral % de P tests 
clasificación binomiales 


B. dia 0,3203 82,2 <0,001 
B. eunomia 0,3089 100 <0,001 











B. selene 0,2931 82,1 <0,001 
B. euphrosyne 0,3566 89,9 <0,001 
B. pales 0,1518 95 <0,001 
B. napaea 0,3011 100 <0,05* 


























Tabla II.— Valores AUC (Area Under the ROC Curve), valor del umbral logístico de la prueba de igualdad de la 
sensibilidad y especificidad de los valores de validación, porcentaje de clasificación y significación de los 11 tests 
binomiales realizados, para cada especie del género Boloria considerada en el estudio. “indica que uno de los once 
tests no fue significativo. 


1.2. Significancia estadística del modelo: Once tests binomiales 

De los 66 tests binomiales realizados para todas las especies del género Boloria, todos excepto los 
de B. napaea fueron significativos con un nivel de p<0,001, con lo que los modelos parecen ser fiables. 
Para B. napaea diez tests presentaron un nivel de significación de p<0,05 y uno de los once no resultó 
significativo. 


1.3. Porcentaje de clasificación 
Se comprobó que más del 80% de las cuadrículas de validación presentaban valores logísticos de 
probabilidad mayores que los requeridos para cada modelo (ver Tabla III). 


Discusión 


La distribución de las mariposas diurnas se ve afectada directa o indirectamente por diversos fac- 
tores, entre ellos la temperatura, precipitación, latitud, altitud (MARTÍN & GURREA, 1990; HAW- 
KINS et al., 2003), o humanos (STEFANESCU et al., 2004). Al analizar la influencia de las variables 
climáticas sobre la distribución de mariposas de montaña, donde los cambios afectan de manera más 
pronunciada (BENISTON et al., 1997), las variables que resultaron más importantes para la creación de 


SHILAP Revta. lepid., 41 (162), junio 2013 273 


H. ROMO, P. SANABRIA & E. GARCÍA-BARROS 


las predicciones en el género Boloria fueron la precipitación del trimestre más cálido y la temperatura 
media del trimestre más seco. Ambas variables están relacionadas con el calor y la disponibilidad de 
agua en verano, momento en que estas especies se encuentran en la fase de adulto. En general, serían 
favorables áreas con temperaturas medias inferiores a 16°C aproximadamente, temperaturas medias 
anuales inferiores a 8°-10°C y precipitaciones en el trimestre más cálido superiores a 150 mm. La ex- 
cepción la constituye B. napaea que prefiere condiciones más frías y húmedas en verano (inferiores a 
10°C y a partir de 250 mm). B. selene es menos selectiva, con mayores probabilidades de presencia a 
temperaturas veraniegas entre 15° a 20°C y precipitación media anual entre 600 y 1.400 mm. De este 
modo se puede observar que con nuestros resultados, las distribuciones potenciales y conocidas del gé- 
nero Boloria se asemejan notablemente, concentrándose en las regiones montañosas de la zona septen- 
trional ibérica, especialmente en la cordillera Cantábrica y los Pirineos. Asimismo, estas zonas coinci- 
den con áreas de mayor diversidad de especies de mariposas (ROMO et al., 2007), zonas menos 
perturbadas y degradadas con gran diversidad de hábitats (LOBO et al., 2001) y donde se han realizado 
los mayores esfuerzos de muestreo (ROMO & GARCÍA-BARROS, 2005). 

A pesar de que se eligieron los valores umbrales para cada especie según el mejor ajuste de las 
predicciones a las presencias ya conocidas, las predicciones obtenidas con el mayor valor umbral no se 
ajustan exactamente a la distribución conocida actualmente en todos sus puntos, incluso algunos se en- 
cuentran en áreas ampliamente desfavorables. Estos puntos se consideran datos atípicos desde el punto 
de vista climático y pueden ser debidos a diferentes causas: errores de georreferenciación en la base de 
datos, identificaciones incorrectas, o debidos a la existencia de otras variables de origen no climático 
que no se están teniendo en cuenta y que pueden tener influencia en la distribución de estas especies. 
En el caso de B. dia, B. selene, B. euphrosyne y B. pales aparecen unas cuadrículas totalmente desfavo- 
rables para la distribución de estas especies y sin embargo proceden de registros en la base de datos. 
Estas cuadrículas se comprobaron a posteriori resultando la mayoría de ellas errores en la base de datos. 
Alguna de las cuadrículas era incorrecta debido a que la altitud de la cita era muy inferior a la de la cua- 
drícula asignada, había sitios mal georreferenciados, posibles errores de identificación de especies de 
alta montaña citadas a 50 m con lo que probablemente no se encontrase allí, o citas muy antiguas en las 
que es posible que la especie estuviese en un área climáticamente desfavorable. De este modo, parece 
que los modelos están siendo bastante fiables, señalando errores de la base de datos en áreas desfavora- 
bles para la supervivencia de la especie. 

En relación al cambio climático, desde finales del siglo XIX se observa que la temperatura mun- 
dial aumenta progresiva e inusualmente (JACOBY & D’ARRIGO, 1997). Este aumento incide sobre 
los patrones de temperatura y precipitación, importantes para la distribución de las mariposas. Al extre- 
marse las condiciones climáticas las distribuciones de las especies estudiadas deberían contraerse incre- 
mentándose su altitud media, en concordancia con la evidencia proporcionada por, por ejemplo, WIL- 
SON et al. (2007) o MERRILL et al. (2008). Las distribuciones predichas sugieren un acantonamiento 
en las áreas pirenaicas y cantábricas, de mayor altitud. B. dia restringe su distribución según los dife- 
rentes períodos considerados, pero mantendría todavía reductos ibéricos importantes en el año 2080, en 
los Pirineos y la cordillera Cantábrica. B. selene sigue una tendencia similar, siendo la segunda especie 
con mayor número de cuadrículas favorables en 2080, pero tiene una tendencia a la pérdida de las cua- 
drículas cantábricas en el escenario más drástico, al igual que B. euphrosyne. Las especies B. eunomia, 
B. napaea y B. pales ya no presentan cuadrículas favorables en la cordillera Cantábrica en 2050, y las 
dos últimas no presentan cuadrículas altamente favorables en 2080. 

A nivel europeo este género presenta la misma tendencia que la aquí mostrada (SETTELE et al., 
2008). Sin embargo, la escala y el área de estudio están influyendo en los resultados obtenidos por los 
modelos y aunque B. selene, B. euphrosyne y B. pales muestren unas tendencias muy similares con las 
obtenidas por nosotros, B. eunomia ni siquiera presenta actualmente áreas favorables en la Península 
Ibérica y B. dia presenta un mejor pronóstico a nivel europeo para el año 2050, aunque termina equipa- 
rándose al obtenido a escala ibérica en el año 2080. 

B. napaea es un caso especial. Es una especie considerada rara en el área Ibero-Balear (al igual 
que B. eunomia) (ROMO et al., 2007). Cuenta con un número muy bajo de registros y se encuentra 
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muy localizada entre los 1.800 y los 2.500 m. A pesar de que esto no supone un problema para Maxent 
aparentemente (HERNANDEZ et al., 2006; PHILLIPS et al., 2006), los resultados obtenidos en los 
modelos son los menos fiables y habría que considerarlos con precaución. Si bien es cierto que el por- 
centaje de clasificación obtenido es el 100%, el 25% de los datos que se usan para validar el modelo 
partiendo de ocho datos de presencia es un número muy pequeño (dos), y es relativamente fácil que 
ambos puntos se hayan clasificado correctamente. Por otro lado, el número de cuadrículas favorables 
para la distribución de esta especie actualmente ya ha sido modelado con anterioridad con diferentes 
programas: Domain, estimadores no paramétricos como ACE, ICE, Chao 1 y 2, Jackknife 1 y 2, Boots- 
trap o un estimador asintótico como es la función de Clench (ROMO et al., 2006) y GARP (ROMO, 
2007). En este caso Maxent propone un mayor número de cuadrículas con condiciones climáticamente 
favorables para la especie, con entre 5 y 23 cuadriculas de diferencia en la predicción, pero siempre 
menor que el número de cuadrículas obtenido con la predicción por GARP (ROMO, 2007). Realmente 
la diferencia obtenida con Maxent con los primeros modelos no es tan grande, ya que en este caso se 
han contabilizado todas las cuadrículas de la predicción, pero es verdad que de ellas, aproximadamente 
12 cuadrículas se encuentran en territorio francés, con lo que los valores se ajustan más a los obtenidos 
anteriormente. En definitiva, habría que intentar favorecer el estudio de esta especie, ya que es una es- 
pecie que necesita unas condiciones muy específicas para su supervivencia siendo altamente vulnerable 
a cambios en su hábitat, y a que no se encuentra en ninguna cuadrícula UTM con más del 15% de su 
superficie dentro de la red de espacios naturales protegidos (ROMO et al., 2007). 

B. dia y B. selene fueron las especies que presentaron un mejor pronóstico en base al cambio cli- 
mático, probablemente debido a que ambas presentan una gran diversidad de hábitats (alta montaña o 
claros o prados de bosque), mientras que las especies de alta montaña (B. pales y B. napaea) presenta- 
ron un panorama dramático, llegando prácticamente a extinguirse sus áreas favorables en 2050 en am- 
bos escenarios. 

Desde el punto de vista de la conservación es una situación inquietante. Se predice que hasta el 
año 2080 dos especies desaparecerán del área de estudio, así como una drástica reducción del área de 
distribución de las demás, lo que concuerda con las predicciones obtenidas a nivel europeo de una gran 
mayoría de especies de mariposas (SETTELE et al., 2008). Cabe pensar si esta disminución de áreas 
favorables predicha por los modelos es suficiente como criterio para introducir a estas especies en listas 
rojas o debería plantearse un plan de seguimiento de las especies. La existencia de áreas climáticamente 
favorables por sí sola no garantiza el futuro de las especies, habría que incorporar factores locales, bió- 
ticos y dinámicos que afecten la distribución de las poblaciones. Pero si además se da el caso de desa- 
parición de esas áreas favorables, la supervivencia de la especie se encuentra aún más amenazada, con 
lo que debería ponerse en marcha un plan de actuación para aplicar medidas de seguimiento de estas es- 
pecies particularmente sensibles, y comprobar en un plazo más corto (10-15 años) si estas tendencias se 
confirman y, en caso afirmativo, proceder a la ejecución de procedimientos más drásticos a fin de pre- 
venir la desaparición de estas especies. 

A pesar de que existe una discrepancia entre diferentes autores sobre si el aumento de temperatura 
debido al cambio climático será más intenso (BENISTON et al., 1997) o no (VUILLE & BRADLEY, 
2000) en zonas de alta montaña, más estudios deben llevarse a cabo para contar con más argumentos y 
poder decidir al respecto (NOGUÉS-BRAVO et al., 2007). Nuestro caso parece apoyar la teoría de que 
las especies de alta montaña van a sufrir más los impactos del cambio climático, mediante una disminu- 
ción en el número de áreas favorables. El aumento en altitud de las especies reducirá cada vez más las 
posibles zonas de colonización, perdiendo las áreas más bajas de su distribución, hasta llegar a agotar 
las regiones de mayor altitud a las que acudir (WILSON et al., 2007). Así pues, y de acuerdo con los 
resultados obtenidos, este factor perjudicará en mayor medida a las especies de alta montaña (alpinas o 
subalpinas), como B. eunomia, B. pales y B. napaea, que se verán afectadas por la falta de hábitats fa- 
vorables en menor tiempo, y a su vez, las especies B. dia, B. selene y B. euphrosyne, que pueden perder 
además los claros de bosque situados a menor altitud de los que disfrutan en la actualidad. 

Los modelos utilizados pueden considerarse válidos con las consideraciones comentadas anterior- 
mente. Es verdad que la decisión sobre cuál es la mejor medida de precisión de los modelos está actual- 
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mente en debate, pero aunque todavía no existe un consenso, el valor del AUC (Area Under a Receiver 
Operating Characteristic (ROC) Curve) es uno de los candidatos con gran aceptación (MANEL et al., 
2001; MC PHERSON et al., 2004; BENITO DEL PANDO & PEÑAS DE GILES, 2007). Véase sin 
embargo LOBO ef al. (2008). La facilidad de interpretación de AUC y su independencia del umbral 
que se considere para identificar la presencia o ausencia de la especie en la distribución, empujan a su 
elección como medida de fiabilidad. Los tests binomiales fueron significativos en la mayoría de los ca- 
sos indicando una buena significación estadística de los modelos. Por otra parte, los porcentajes de cla- 
sificación que se obtuvieron fueron superiores al 82% en todos los casos, por tanto una gran cantidad 
de las cuadrículas de validación mostraban valores correctos de predicción, recalcando la fiabilidad de 
los modelos. 

Con el fin de aumentar la posibilidad de predecir y prevenir los efectos del cambio climático sobre 
la biosfera a corto o largo plazo, será necesario continuar realizando estudios que permitan conocer ca- 
da vez más y mejor lo que está sucediendo en realidad. En este caso, sería de gran ayuda ampliar el es- 
tudio utilizando otras variables no consideradas en el presente trabajo como por ejemplo las distribucio- 
nes de sus plantas nutricias, variables históricas u otras variables ambientales que puedan tener 
influencia sobre su distribución. Además, debido a su distribución restringida y a su susceptibilidad al 
cambio climático, se hace necesario el seguimiento de las poblaciones del género Boloria con el fin de 
comprobar que no cambie su estatus y necesiten ser incluidas en alguna lista roja. 
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Figura 1.- Distribución potencial actual realizada con el programa MaxEnt (formato logístico) para las espe- 
cies del género Boloria en la península Ibérica. A: B. dia; B: B. selene; C: B. euphrosyne; D: B. eunomia E: B. 
pales; F: B. napaea. Los tonos negros indican áreas adecuadas ambientalmente para la presencia de cada espe- 
cie, disminuyendo esta favorabilidad desde los tonos de gris más oscuros a los más claros. Los puntos blancos 
representan la distribución actual conocida de las especies. 
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Figura 2.— Modelo de distribución potencial obtenido con el programa MaxEnt (formato logístico) proyectado 
bajo el escenario HadCM3 de cambio climático para el año 2020 (arriba), 2050 (centro) y 2080 (abajo), para la 
especie Boloria dia en la península Ibérica. El valor de umbral más alto donde se encuentran las cuadrículas cli- 
máticamente más favorables se representa en negro, disminuyendo su valor según la escala de grises. 
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Figura 3.— Modelo de distribución potencial obtenido con el programa MaxEnt (formato logístico) proyectado 
bajo el escenario HadCM3 de cambio climático para el año 2020 (arriba), 2050 (centro) y 2080 (abajo), para la 
especie Boloria selene en la península Ibérica. El valor de umbral más alto donde se encuentran las cuadrículas 
climáticamente más favorables se representa en negro, disminuyendo su valor según la escala de grises. 
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Figura 4.— Modelo de distribución potencial obtenido con el programa MaxEnt (formato logístico) proyectado 
bajo el escenario HadCM3 de cambio climático para el año 2020 (arriba), 2050 (centro) y 2080 (abajo), para la 
especie Boloria euphrosyne en la península Ibérica. El valor de umbral más alto donde se encuentran las cuadri- 
culas climáticamente más favorables se representa en negro, disminuyendo su valor según la escala de grises. 
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Figura 5.— Modelo de distribución potencial obtenido con el programa MaxEnt (formato logístico) proyectado 
bajo el escenario HadCM3 de cambio climático para el año 2020 (arriba), 2050 (centro) y 2080 (abajo), para la 
especie Boloria eunomia en la península Ibérica. El valor de umbral más alto donde se encuentran las cuadrícu- 
las climáticamente más favorables se representa en negro, disminuyendo su valor según la escala de grises. 
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Figura 6.— Modelo de distribución potencial obtenido con el programa MaxEnt (formato logístico) proyectado 
bajo el escenario HadCM3 de cambio climático para el año 2020 (arriba), 2050 (centro) y 2080 (abajo), para la 
especie Boloria pales en la península Ibérica. El valor de umbral más alto donde se encuentran las cuadrículas 
climáticamente más favorables se representa en negro, disminuyendo su valor según la escala de grises. 
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Figura 7.— Modelo de distribución potencial obtenido con el programa MaxEnt (formato logístico) proyectado 
bajo el escenario HadCM3 de cambio climático para el año 2020 (arriba), 2050 (centro) y 2080 (abajo), para la 
especie Boloria napaea en la península Ibérica. El valor de umbral más alto donde se encuentran las cuadrícu- 
las climáticamente más favorables se representa en negro, disminuyendo su valor según la escala de grises. 
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